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We investigate the effect on biomembrane mechanical properties due to the presence an external 
potential for a non-conductive non-compressible membrane surrounded by different electrolytes. 
By solving the Debye-Hiickel and Laplace equations for the electrostatic potential and using the 
relevant stress-tensor we find: in (1.) the small screening length limit, where the Debye screening 
length is smaller than the distance between the electrodes, the screening certifies that all electrostatic 
interactions are short-range and the major effect of the applied potential is to decrease the membrane 
tension and increase the bending rigidity; explicit expressions for electrostatic contribution to the 
tension and bending rigidity are derived as a function of the applied potential, the Debye screening 
lengths and the dielectric constants of the membrane and the solvents. For sufficiently large voltages 
the negative contribution to the tension is expected to cause a membrane stretching instability. For 
(2.) the dielectric limit, i.e. no salt (and small wavevectors compared to the distance between the 
electrodes), when the dielectric constant on the two sides are different the applied potential induces 
an effective (unscreened) membrane charge density, whose long-range interaction is expected to lead 
to a membrane undulation instability. 



I. INTRODUCTION 



Bio-membranes are thin fluid films composed mostly 
of lipids. In cells they help separating different cellu- 
lar environments and compartments. Biomembranes are 
typically "soft", i.e., the typical energy required to bend 
them is of the order the thermal energy and membrane 
tension is often quite small. Softness implies that mem- 
brane geometry can become sensitive to different per- 
turbations, such as alteration of the electrostatic con- 
figuration. Much effort has for instance been devoted 
to calculation of the electrostatic contribution to tension 
and bending rigidity for membranes with fixed charges or 
surface potentials in an electrolyte solution, see [l[ for a 
review; in general the presence of fixed (and screened due 
to the electrolyte) charges tend to increase bending rigid- 
ity and hence make the membrane stiffer. However, it has 
also been found that when the membrane charges are not 
fixed but free to rearrange themselves on the surface and 
no electrolyte solution is present to screen the interac- 
tion, a long-wavelength undulation instability can occur 
0, somewhat similar to DNA condensation Elec- 
tric fields can also be present across intrinsically neutral 
membranes. An example would be a nerve cell, where 
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ion pumps create a potential difference between the two 
sides of the nerve cell membrane Q • Another example is 
provided in laboratories by the routine formation of lipo- 
somes in a process known as electroformation [5J , during 
which lipid membranes are swelling from electrodes un- 
der the application of electric fields. 

In this study we investigate the electromechanical cou- 
pling of a membrane and an applied potential. In particu- 
lar, we solve the Debye-Hiickel and Laplace equations for 
the electrostatic potential for a non-conductive, incom- 
pressible membrane between two fiat electrodes (kept at 
fixed potentials). On either side the membrane is sur- 
rounded by different electrolyte solutions. From the so- 
lutions for the potential we quantify how the correspond- 
ing induced membrane charges change the free energy for 
the membrane and identify electrostatic contributions to 
membrane mechanical parameters. In the presence of an 
electrolyte (in the small screening length limit, see be- 
low) we find that the electrostatic contribution to mem- 
brane bending rigidity is positive. In the absence of 
added salt (the dielectric limit) the membrane becomes 
unstable against long wavelength undulations (in a some- 
what similar fashion to the behaviour of the interface be- 
tween two immiscible fluids, see Ref. @) if the two fluids 
surrounding the membrane have different dielectric con- 
stants. For the symmetric dielectric well as for 
the small screening length limit, the membrane tension 
receives a negative contribution; for a sufficiently large 
applied potential this contribution would lead to a mem- 
brane stretching instability. 

Several studies (see [l| and references therein) have in- 
vestigated the electrostatic contributions to tension and 
bending rigidity for membranes having a fixed surface 
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charge density or fixed potential at the membrane. How- 
ever, less work has been dedicated to the effect of induced 
charges due to an applied potential, as considered in the 
present study. The results found here complement previ- 
ous results given in [7| (using coupled hydrodynamical- 
electric field equations, similar to [8()) for the symmetric 
case of a membrane surrounded by identical electrolytes 
(in the small screening length limit), to the asymmetric 
case and by giving an explicit expression for the bend- 
ing rigidity; in the limit of identical electrolytes on the 
two sides of the membrane our expression for the ten- 
sion agrees with that derived in 0. Also, our formula- 
tion allows us to investigate the dielectric (no salt) limit, 
which was not considered in Q. However, unlike 0] we 
do not consider any dynamical effects. In contrast to 
the study [9( based on electrolyte conductivities, our ap- 
proach through the Debye-Huckel equation allow us to 
study the effect of a finite Debye screening length. In 
Ref. [9( a non-zero membrane conductivity was consid- 
ered, whereas we consider non-conductive membranes. 
Similar to the results in Q we find a negative contribu- 
tion to the membrane tension in the presence of an elec- 
trolyte, although the two results are difficult to compare 
because of the different mathematical formulations. 

This work is organized as follows: In Sec. [ill we give 
the general equations governing the electrostatic response 
of a non-conductive membrane of any shape; the mem- 
brane region is described by the Laplace equation and 
the electrolyte solutions on either side satisfy the Debye- 
Hiickel equations. In a standard fashion, the boundary 
conditions are that the potential and the displacement 
fields should be continuous. In Sec. |Hl] these equations 
are solved for the case of a flat membrane in an exter- 
nal potential. In Sec. IIVI corrections to the flat case 
solutions are derived for a weakly curved incompressible 
membrane. In Sec. [V] the forces acting on the mem- 
brane, as well as the corresponding electrostatic contri- 
bution to the membrane free energy, are obtained. As- 
suming that membrane fluctuations occur on a time scale 
slower than the relaxation time for the electrostatic po- 
tential we use the expressions for the forces for the weakly 
curved membrane in order to obtain the renormalized 
membrane mechanical parameters, such as tension and 
bending rigidity, (in terms of a power series expansion in 
the wavevector) as a function of the applied potential, the 
salt concentrations (entering through the Debye screen- 
ing lengths) and the dielectric constants of the membrane 
and the solvents. We investigate three different limits: 
(1.) the small screening length limit, where the Debye 
screening length is smaller than the distance between the 
electrodes; (2.) the dielectric limit, i.e., no salt; (3.) the 
symmetric case, where the salt concentration and dielec- 
tric constants on the two sides of membrane are equal. 
The results for the membrane mechanical parameters in 
the three limits above (the main results in this study) are 
given in Eqs. (j¥T|) - (|5ip . Finally, in Sec. I VII a summary 
and discussion are given. 



II. GENERAL FORMULATION 

We are interested in how biomembrane mechanical pa- 
rameters (and thereby, for instance, membrane fluctu- 
ations) are effected by an applied potential. Two pa- 
rameters characterizing a membrane in the absence of 
an applied pot ential are the tension a and the bending 
stiffness K jlOl . [Tl| . As an external potential is applied 
there will in general be electrostatic contributions (<7 e i 
and K c \) to both of these quantities so that a — > a + a c \ 
and K —* K + K c \ in the presence of the applied poten- 
tial. An aim of this paper is to calculate a e \ and K c \. 
In a standard fashion we consider a small perturbation 
from a flat membrane, characterized by a height undula- 
tion h(x, y) (where x and y are coordinates in the plane 
of the flat membrane) and solve the electrostatic equa- 
tions (via a Fourier-transformation in the x- and y coor- 
dinates) for this weakly perturbed geometry. Through a 

power series expansion in wavevector q (q = ^Jq 2 + q? , 

where q x and q y are the Fourier-transform variables of x 
and y respectively) of the free energy G one may iden- 
tify cr e i and K e \ (see chapter 2 in Ref. (HI). We find it 
convenient to, rather than utilize the free energy directly, 
consider the electrostatic contribution to the "restoring" 
force, from which we identify a c \ (K c i) as the prefactor in 
front of the —q 2 h (—q h) term in a small g-expansion of 
the force, where h = h(q x , q v ) is the Fourier-transform of 
h{x,y). Notice that for obtaining the tension and bend- 
ing rigidity it suffices to keep terms linear in h in the 
restoring force expression. In general there may be other 
terms in the power series expansion in q. For instance, as 
noted in the introduction, for the asymmetric dielectric 
case there is a membrane undulation instability which 
mathematically arises due to the presence of a negative 
term linear in q in the series expansion. 

The approach described above relies on a "quasi- 
static" approximation, i.e. we assume that membrane 
fluctuations occur on a time scale t mcm slower than the 
time scale t \ over which the electrostatic configuration 
adjusts itself (t e \ <C t mcm ). This assumption allows us to 
solve the electrostatic problem for a fixed, weakly curved 
(but otherwise arbitrary) geometry. Let us estimate the 
time scales t c \ and t mcm : To estimate t c \ for an elec- 
trolyte we assume that this time scale equals the time 
for an ion to diffuse the distance of the order the Debye 
screening length, i.e. t c \ as k~ 2 /D, where D is the ion 
diffusion constant and k is the inverse Debye screening 
length introduced below (one may more realistically as- 
sume that t \ = min{K~ 2 , q~ 2 }/D, since for a wavelength 
perturbation of the order 1/q the ions need only diffuse 
a distance 1 jq for the ion cloud to relax) . From the Ein- 
stein relation and Stoke's law, we have D = ksT/ (6iri]R), 
where rj is the viscosity, fcs is the Boltzmann constant, 
T the temperature and R the ion Stoke's radius, using 
k B T = 4 • 1CT 21 J, rj = 1(T 3 Ns/m 2 and R « 0.1 - 0.3 
nm [l2j], we find D w 10~ 9 m 2 /s. Furthermore, taking 
= 10 nm we obtain: t c \ « 10~ 7 s. For the membrane 
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relaxation time we estimate 13 1 t mcm ~ r]/(q 3 K), and as- 
suming q- 1 > 100 nm, K = 10~ 19 J we find t mcm > 10~ 5 
s. Therefore, indeed, we have t c \ <C t mcm in general. 
We note from the expressions above that the longer the 
wavelength perturbation (smaller q) the better justified 
is our quasi-static approximation. In the dielectric limit 
there are no ions and the relevant relaxation time t e \ is in- 
stead that of water relaxation (hydrogen-bond rearrange- 
ment time), which typically is of the order 10" 12 s (see 
Ref. [HI) at room temperature, again certifying that 

^cl ^mem* 

We are now set to consider the the effect of an applied 
potential on the mechanical properties of a biomembrane 
within the quasi-static approximation. The explicit prob- 
lem we consider is depicted in Fig. [T] an incompressible 
membrane of thickness 2d is placed with its center-of- 
mass at positions z = between two flat electrodes (at 
z = ±(L + d)) which are kept at potentials =fA(/>/2; the 
distance between the membrane surface (for a flat mem- 
brane) and the electrodes are hence L. Regions 1 and 
3 are electrolyte solutions, and in general these two re- 
gions are of different composition (different concentration 
of ions and different dielectric constants) . As noted above 
the first stage towards calculating electrostatic forces on 
the membrane and thereby the membrane free energy in 
the presence of the applied potential is to obtain the elec- 
trostatic potential <&(af). In this section we give general 
equations determining <fr(x) for any membrane shape. In 
the subsequent sections we analyze in detail: (i) the flat 
membrane case, see FigOJi); all quantities for this case 
carry a superscript (0); (ii) for the weakly curved situa- 
tion (since we assume the membrane to be incompressible 
we only consider undulation deformation modes), illus- 
trated in Fig. [TJd), there will be corrections of all quan- 
tities compared with the flat case; all such corrections 
carry the superscript (1). 

For all three regions the electrostatic potential satisfies 
the Poisson equation (using SI units [151 ]) 



e e 7 V 2 <i> 7 (f) = -p 7 (x). 



(1) 



We will henceforth use a subscript 7 (=1,2 or 3) to dis- 
tinguish quantities in the three different regions. Above 
£0 is the permittivity of vacuum, e 7 is the dielectric con- 
stant for region 7, and p 7 (x) is the free charge density 
(the bound charges are, in a standard fashion, taken into 
account through e 7 ). We require 



Pj(x)d 3 x = 0, 



(2) 



i.e, that in each region we have charge neutrality (the 
membrane is assumed to be impermeable to ions). 

Let us now consider the explicit expression for the 
charge density in each of the three regions. 
Region 2. In region 2 we assume 



P2(x) = 0, 



(3) 



i.e. there are no free charges in this region. 

Region 1. In region 1 we assume that there are 
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FIG. 1: Cartoon of the problems considered in this study: a 
membrane of width 2d is placed between two electrodes. Re- 
gion 1 and 3 are characterized by dielectric constant £1 and 
£3 and Debye screening lengths Ki and respectively. The 
membrane is assumed non-conductive (and non-compressible) 
and characterized by a dielectric constant £2. Solving for the 
electrostatic potential $(x) for (a) a flat membrane, and (b) 
a weakly curved membrane (and utilizing the stress-tensor) 
allows us to obtain the electrostatic contribution to the mem- 
brane mechanical parameters. 



ions of two types (positively and negatively charged) 
which are taken to be Boltzmann-distributed, i.e. 
(i distinguishes the different ionic species) pi{x) — 
J2i 9 Ic buik,i exp[-/3g*($i(^) - <f>f)] where q l is the charge 
of ionic species i, c^ ulk j is the concentration of ions in 
region 1, and f3 = 1/QcbT), with ks the Boltzmann 
constant and T the temperature, as before. The con- 
stant (pi is determined through the charge neutrality 
condition. Making a linear approximation, i.e. assum- 
ing (3q l ($i(x) — <Pi) <§C 1, and using the charge neutral- 
ity condition in the absence of the external potential, 
Si 9* c buik = 0, we find that the charge density can be 
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written: 



with 



pi(x) = -e £i«i(*i(s) - 0?), 



£q£i 



(4) 



(5) 



where ki is the inverse Debye screening length for region 
1. 

Region 3. For region 3, with the same approximations as 
above, the charge density becomes: 



where 



p 3 (x) = -e e 3 Kl(<f> 3 (x) - 4>° 3 ), 



K 3 — ~~ y~l[g'] 2c bulk.3 

£0^3 ; 



(6) 



(7) 



where h is the normal to the respective interface. Eqs. 
(O, (JHJ) and (HHJ), together with the boundary conditions 
Eqs. (jlip. (|12p and (fl~3f completely determine the elec- 
trostatic potential 

From the solutions for $(x) one can calculate other 
quantities. For instance, one may obtain the induced 
potential defined through < &i n d(a?) = — && PP \{x) , 

where the applied potential is & app i(x) = —A<fiz/[2(L + 
d)]. The total electric field is given by E(x) — — V$(x), 
the applied electric field is E app i = — V$ app i and the 
induced field is E m d = — V$i n d = E — E app \. In the next 
section we solve the equations given in this section for 
the case of a flat membrane. In the section after we find 
corrections to the flat case solutions for a weakly curved 
membrane. 



III. POTENTIAL FOR A FLAT MEMBRANE 



is the square of the inverse Debye screening length and 
c k>uik 3 i s the concentration of ions in region 3. 

Inserting Eqs. {3]) into Eq. (JTJ yields the Laplace equa- 
tion 



V 2 $ 2 (f) = 



(8) 



for the potential in region 2. Inserting (HJ) and into 
Eq. ([1]) gives the Debye-Hiickel equation (7 = 1,3) 



V 2 $ 7 (x) - k*($ 7 (2) - 4>X) = 



(9) 



for regions 1 and 3. The charge neutrality condition, Eq. 
@, for region 1 and 3 becomes 



[$ 7 (f) -c^d^x = 0. 



(10) 



Note that we have trivial charge neutrality in region 2 
[seeEq. ©]. 

Let us now consider the boundary conditions supple- 
menting the equations above. At the electrodes we have: 



®l(x)\g=-L-d 
$3(x)\z=L+d 



A0 

A<f> 
2~' 



(11) 



In addition we have that the potential and the nor- 
mal component of the displacement fields are continuous 
across the reg ion 1-region 2 and region 2-region 3 bound- 
aries, i.e. [l5( 

<J>i(af) = $2(x) & t region 1 — 2 boundary, 
^2(2;) = ®3(x) at region 2 — 3 boundary (12) 

and 

£ifi ■ V$i(x) = eifi ■ V$2(#) at region 1 — 2 boundary, 
£2n ■ V$2(a;) = £3" • V$3(x) at region 2 — 3 boundary, 

(13) 



Below we obtain the electrostatic potential in and 
around a flat membrane. We use a superscript (0) to 
indicate the flat case quantities. 

For a flat membrane the solutions depend only on z, 
and explicitly the solutions to Eqs. ijHJl and © are 



$?\z)=cp 2 +A 2 z 



(14) 



where (j)® and A2 are constants determined by the bound- 
ary conditions below. Also, 

<pf\z) = 0° + A x {e^ s+d) - e - Klt - L+d+z ^ (15) 

and 

<4 0) (z) = 4>l + A 3 {e~ K3{z - d) - e -K3( L +d-*)) (16) 

where (j>°, A\ and A 3 are constants, and we used the 
charge neutrality condition, Eq. (|10[) . 

We now use the boundary conditions (together with 
the fact that the boundary surfaces are at z = ±d for 
the flat case considered here, see FigQJ,) in order to de- 
termine the unknown constants above. From Eq. (jlip . 
the condition that the potential is continuous, Eq. (TT2"]) 
and the fact that the displacement field is continuous, 
Eq. (fT3"|) , we get 6 equations for the 6 constants c/> 7 and 
A 1 (7 = 1, 2, 3). Solving these equations leads to 

A(f> 



A, = 



2(1 + e- 

2e 2 
A<f> 



-K\L 



2(1 + e- K * L ) 



l 3 T 



(17) 



and 



—^(g{Ki)h -g(K 3 )k)r 
-^(i-.W) 



(18) 
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where 



/(<?) = 



1 - e- qL 



1 



= tanh( 



(19) 



and r = \/[g{ni)li + .9(^3)^3 +d/e 2 ], and we introduced 
the "rescaled" Debye screening lengths li — (eiKi) -1 and 

^3 = (£3«3 

There are three limits of particular interest: 

1. "small" screening length, kiL,k^L 3> 1. For this 
case we have g(q) — » 1 in Eqs. (JTTJ) and (TTSJ). Also 
the prefactors for A\ and A3 simplify. 

2. we define i/ie dielectric limit as the limit of no salt, 
i- e - c buik i' c buik 3 ~* 0; within the Debye-Huckel 
approximation this is the equivalent to m, K3 — > 0. 
Expanding the exponentials in Eqs. (TT5)) and (HHJ), 
and using the explicit form for A\ and A3 above 
one straightforwardly show that the solutions in all 
three regions take the form <f>(x) = az + b where a 
and b are constants independent of ki and K3, as 
it should since in the dielectric limit the potential 
satisfies Laplace equation, V 2 $(x) = 0, see Eq. ((9|. 



3. the symmetric case, 
limit we find that A\ 



£\ = £3 and Ki 
= -A 3 , 4>\ = -< 



= k 3 . In this 
>3> and <t>l = 0. 



The potential and charge densities are illustrated in 
Fig. [2j using the flat membrane results in Eqs. (fT4]) , 
(fT"5|) and l|16|> . The electric field in the z-direction (the 
electric field components in the x- and y-direction are 
zero for a flat membrane) is also illustrated. We notice 
that the potential is continuous as it should and that 
the free charges tend to build up close to the membrane 
and electrodes (for ki,k 3 ^ 0). Since the normal com- 
ponent of the displacement field is continuous across the 
boundaries, the relative jump in the electric field as the 
the boundary between region 1-2 (region 2-3) is crossed 
equals £2/^1 (£2/53), see Fig. [5] (bottom). 




-0.5 



S2 





z/L 



0.5 




— E 

E 



appl 




-0.5 





z/L 



0.5 



IV. POTENTIAL FOR A WEAKLY CURVED 
MEMBRANE 

We now consider a weakly curved membrane, see Fig 
[]J>): the center of the membrane is slightly displaced 
from the flat (z = 0) case, according to z = h(x, y) with 
membrane surfaces at z = ±tf + h(x,y) (l6j . We write 
the solution for the electrostatic potential according to 
(7 = 1, 2, 3 as before): 

%(x)=$W(z) + $W(x), (20) 

where $7° (z) is the potential for region 7 for the flat 
case given in the previous section and $iy (2?) is a cor- 
rection to the potential due to the perturbed geometry. 
In this section we will calculate only to first order in 
the perturbation h(x,y); this suffices for obtaining mem- 
brane mechanical parameters such as tension and bend- 
ing rigidity (see discussion at the begininning of Sec. HI]) . 



FIG. 2: Illustration of the response of a flat membrane, sur- 
rounded by different electrolyte solutions, to an applied po- 
tential, a) Electrostatic potential <&(af) = & °' (z) for a flat 
membrane (solid curve) as a function of position z. Illustrated 
are also the applied potential $ apP i(z) (dotted curve) and the 
induced potential $i n d(^) (dashed curve); b) Free charge den- 
sity p(z) (solid curve) as a function of position z. The dashed 
curve at p = corresponds to the free charge density in the 
absence of the applied potential (A<f> = 0); c) The electric 
field E = —d§/dz in the z-direction for different positions 
z. Illustrated are also the applied field -B app i = — 9<l?a PP i/9z 
and the induced field -Emd = — d&i n <\/dz. The vertical lines in 
each figure corresponds to the membrane surfaces. The cho- 
sen parameter values are the same for all three graphs, and 
are listed in the top graph (we chose the somewhat unrealisti- 
cally large value £2 = 15 in order to be able to better visualize 
the Debye screening layers, a more reasonable value would be 
£2 ~ 2). Note that the potential is continuous, whereas the 
free charge density and the z-component of the electric field 
are not continuous, as it should. 
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In each of the regions Eqs. (JH) and ([9]) has to be satisfied. 
$ 7 (x) satisfies these equations for h(x,y) = 0, and we 
therefore require that the corrections satisfy: 



V 2 4 1} (f) = 



and (7 = 1, 3) 



V 2 $«(x)- K 2 $«(x) = 0. 



(21) 
(22) 



Let us now consider the boundary conditions (the bound- 
ary surfaces are at z = ±d + h(x, y) for the curved mem- 
brane considered here). Since $ 7 0) (z) satisfies Eq. dJ| 
we require for the perturbation: 



$P(x)U = _£_ d = 0, 



L+d 



= 0. 



(23) 



Before imposing the conditions that the potential and the 
displacement fields are continuous, Eqs. (fT2"j) and (fTB")) . 
we note that for "small" h any scalar quantity g{x) may 
be expanded, to first order in h, according to: 

g(x) = g i0) (x)\z=±d+h + g {1) (x)\z=±d+h 

we briefly discuss the quantitative meaning of "small" h 
at the end of this section. Eqs. (fl~2| and (fT3| can then 
be written 



<9 " 



and 



d 2 ^ } d$\ 

dz 2 
d 2 ^ ] 



e 2 (h 



(i) 



dz 



+ 



dz > 

dz 



dz 

z=-d = £-2{h 
)\ z= d = e 3 {h 



a 2 $ 2 0) 




dz 2 


h 

dz 


d 2 ^ 




dz 2 


1- 

dz 



\z= — di 



and (7 = 1, 3) 



(i) 



(i) 



0, 



where g = J q 2 + q 2 and 



(28) 



(29) 



The boundary conditions Eqs. ([23|). ([24| and d25j) re- 
mains the same in Fourier-space (using the fact that $ 
is independent of x and y), with the sole replacement 



(0) 



W _> $W an d ft -> ft. The solution of Eq. (JSTJ) is 

#i 1) (g,^)=C 2 ( 9 )e^ +i? 2 (g)e-«*, (30) 



with g-dependent coefficients C 2 (g) and -D 2 (g) ■ Using the 
boundary condition Eq. (|23[) we find that the solutions 
to Eq. (J28]) are: 



^ 1} (g,z) = £>i( g )(e 9_1 ^ - e -^+ d+2L )) (31) 



and 



$i 1J ( 9 ,z) = £> 3 (g)(e- 



-q 3 (z-d) _ e q3(z-d-2L)\ 



(32) 



The unknown coefficients C2(q), Di(q), D2(q) and D 3 (g) 
are determined through the boundary conditions Eqs. 
H and ©. We find 



D 2 {q) = -h 

e qd k+(p 1 + siq x r{qi)rrq) - e"« d fcf (p 3 - e 3 q 3 r(q 3 )m+) 



e 2qd k^k^ - e- 2 i d k^k~ 3 



(33) 



(25) 



where we introduced the short-hand notations (7 = 1,3) 

k± = k±(q) = e 2 q±e 1 q 1 r(q 1 ), 
m± = A 2 ±A 1 K 1 (l + e-^ L ), 
p 7 = e 1 A 1 K 2 {l - e~ fi -< L ), 



where we used the fact that $7 satisfies the boundary 
conditions for h = 0. Eqs. (|2"Tj) and ([22| together with 
the boundary conditions Eqs. (f2"3"| . (J2H) and (|23|) com- 
pletely determine the correction $7 (x). 

We proceed by introducing the Fourier-transform in 
the x- and y-direction (not for ^-direction) of ^^''(x): 



r(g) 



1 + e~ 2qL 

1 _ e -2qL ■ 



(34) 



with Ay given in Eq. (fl~7|) . The remaining coefficients 
are 

e qd 

C 2 (g) = — (D 2 {q)e qd k+ + h Pl +he 1 q 1 r(q l )m^), 



$ J 1 \q x ,q y ,z)= f dxdye iq * x+iq yV$M(£), (26) ^(g) 



^ _ e -2q 1 L 



and similarly we denote by h(q x ,q y ) the Fourier- 
transform of h(x,y). Eqs. (f2"T) and (|2"2"|) can then be 
written: 



D 



(hm^ +C 2 (q)e- qd + D 2 (qy d ), 
(?) = ! _ l-2- q3L {hmt + C 2 (q)e qd + D 2 (q)e- qd ) . 



(27) 



(35) 

We notice that Ci(g), -Di(g), D 2 (q) and D 3 (q) are all pro- 
portional to ft. as they should be. The full solution for the 
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electrostatic potentials 3> 7 (a?) for a weakly curved mem- 
brane is given by Eq. (20]), where <I> 7 0) (z) were 
the previous section and the Fourier-transform of ^>\ 1> (x) 



given m 



are given by Eqs. (]3"0]). ([51]) and (32]) together with Eqs. 
(33j), (05) and (J35]). 

Again, we point out that in order to obtain mem- 
brane mechanical parameters, such as tension and bend- 
ing rigidity, it suffices to know the restoring force (cal- 
culated in the next section) to first order in h, i.e. it 
is enough to consider "small" fluctuation amplitudes. 
Therefore, even thought the results given above formally 
assume that h is smaller than all other length scales in 
the problem (d, K~ , q^ 1 and L), they are sufficient for 
the purpose of calculating the tension and bending rigid- 
ity. A different matter is whether the Helfrich form of the 
electrostatic contribution to the restoring force (in terms 
tension and bending rigidity), derived in the next section, 
describe well "large" membrane fluctuations in the pres- 
ence of an applied potential. To address this question we 
must clarify the meaning of "small" h, i.e. make clear 
what is the relevant dimensionless expansion parameter, 
in the context of calculating the electrostatic contribu- 
tion to the membrane restoring force or membrane free 
energy. The electrostatic contribution to the membrane 
free energy is determined by the interactions between in- 
duced charges at or in the vicinity of the membrane, with 
characteristic interaction distances of the order d and k" 1 
(here we consider the small screening length limit, where 
K-y is non-zero). Therefore, whenever dq <C 1, n^ 1 q *C 1 
and hq <C 1 all interactions are effectively local on a lo- 
cally flat membrane and the free energy must take the 
Helfrich form [l(| • Note that this argument is valid inde- 
pendent of particular values of /ik 7 and h/d, and the free 
energy expansions given in the next section (in the small 
screening length limit) is therefore an expansion in the 
(small) parameters hq, dq, K~ 1 q, but when these parame- 
ters are small there is no restriction on the values of /ik 7 
and h/d. In consistency with the discussion above, we 



point out that in Ref. [17| a method that formally avoids 
the assumptions /ik 7 <C 1 and h/d <C 1, by utilizing the 
geometrical transformation z' — z — h(x,y), was found 
to give results consistent with results obtained by the 
flat membrane perturbative approach (of the kind used 
in this paper) for a charged membrane in an electrolyte. 



V. FORCES, FREE ENERGY AND 
ELECTROSTATIC CONTRIBUTION TO 
MEMBRANE MECHANICAL PARAMETERS 

In this section we derive general expressions for the 
membrane forces (within the quasi-static approximation) 
and the corresponding electrostatic contribution to the 
membrane free energy, using the results from the previous 
two sections. In particular, we obtain the electrostatic 
contribution to the membrane free energy for three cases: 
(1.) the small screening length limit, where the Debye 
screening length is smaller than the distance between the 



electrodes; (2.) the dielectric limit, i.e., no salt, and; (3.) 
the symmetric case, where the salt concentrations and 
dielectric constants on the two sides of the membrane 
are equal. 



A. Membrane forces via a stress-tensor calculation 

The forces acting on the membrane are obtained using 
the relevant stress-tensor, Tj,-. Following the derivation 
in Appendix[X]we have (in the Debye-Hiickel regime con- 
sidered here) 



T — 



fd<t> y <9$ 7 



7 1 dxi dxj 



7 "^7 

dx k 



dx k 



■plS. 



7"u > 



(36) 



where x\ = x, xi = y and X3 = z. The first two terms are 
the usual Maxwell stress-tensor [ID, [l8| , the third term 
is an osmotic contribution for the ions being "confined" 
by the electric potential (see appendix [X]) , and the last 
term incorporates pressures p® for each of the three re- 
gions (7 = 1,2,3). The discontinuities of Ty- at the region 
boundaries will produce forces on the membrane which 
will have to be balanced by other forces in the system, 
such as for example viscous forces within the membrane 
or from the surrounding bulk fluids. We are only inter- 
ested in calculating the electromechanical contribution at 
a given (x, y) to this total force balance here. To do this 
we note that the force (per unit area) in the «-direction 
on a region boundary from the stress in a given region is 
± n jTji evaluated at the boundary, where the mem- 
brane normal rij is taken to point towards positive z and 
the plus (minus) sign applies when the region is at larger 
(smaller) z than the boundary. Defining 



f(0) 
J7 



4 0) 



4 X) 



1 <9$ (0) 



(0) 

7 



\z=±d 



-P7 



,9$ 7 0) a 2 $ 7 0) 

■ [n- 



(i) 



dz 



dz 2 



dz 



(0) 



dz 



(37) 



where z = —d (z = d) is to be used for forces acting on 
the interface separating region 1 and 2 (region 2 and 3), 
and using that to first order in h: n z = 1, rij = —djh 
(j = x,y), we find that the z-component of the total 
force acting on the surface separating regions 1 and 2 is 
fi—2 — /2 — /i- Using the explicit expressions for the 
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potentials from the previous two sections we find 

h-2 = fi% + A% 

f[% = - £0 ( 2ei Af K fe-^-^i)+rf-^ 

f[% = -eo(^i(g)eiA lKl [g 1 (l + e - Kli )(l + e- 2 « li ) 
-/si(l - e~ KlL )(l - e~ 2 ^ L )\ 



-s 2 A 2 q[C 2 (q)e-^ - D 2 (q)e" d } 



(38) 



where /}_2 ^ s the f° rce on a &at membrane interface, 

and fi^2 1S the first order correction for a weakly curved 
membrane, here expressed explicitly in terms of its 
Fourier-transform f[_ 2 ; note that f[_ 2 is proportional 
to h (via D\, C 2 and D 2 ), as it should. Similarly, the z- 
component of the force acting on the surface separating 
regions 2 and 3 is f 2 —z — Jz — f 2 and we find 



/2-3 — 

f (o) 

J 2 -3 ~ 

J(l) _ 

J 2-3 ~ 



AO) , 

./2-3 "T" J 2-3 

eo(2e 3 ^3«3e" Ks 



£ (D 3 (q)e 3 A 3 K 3 [q 3 (l + e- K3L )(l + e" 

-K 3 (l-e- K3L )(l-e- 2 « sL )} 
-s 2 A 2 q[C 2 (q)e« d - D 2 (q)e-* d }) (39) 



where A 7 are given in Eq. (fI7| and C 2 (q), D 1 {q) are 
given in Eqs. (f3"3"]) and ([33)1 . The results given in Eqs 
(|38| and p9|) completes the calculation of the forces act- 
ing on the membrane interfaces. In appendix [Bl wc utilize 
these results in order obtain results for the total net force 
on a flat membrane in some detail. In the next subsection 
the membrane free energy and the electrostatic contribu- 
tion to the membrane mechanical parameters in different 
limits are investigated. 



the free energy and electrostatic contribution to mem- 
brane mechanical parameters (through a power series in 
q, i.e. a long wavelength expansion). In particular we 
want to compare the results of such an expansion to the 
corresponding result for a "free" membrane: the free en- 
ergy G for a membrane is described by the Helfrich form 
[la, [n| G = JdA {2KB 2 + a) where H is the mean 
curvature, dA the area element on the membrane, a is 
the tension and K is the bending rigidity. The restoring 
force is then obtained as f IS = —6G/5h giving in q-space 



f TB (q) = -[aq 2 + Kq i ]h + 0(h 2 ). 



(40) 



This type of expansion requires only that the expectation 
value of (yh(x,y)) 2 is small [ll|, i.e. that the character- 
istic fluctuation amplitude is small compared to l/q. In 
the presence of an applied potential there will be elec- 
trostatic contributions a c \ and K e \ to the tension and 
bending rigidity, so that a — > a + a e i and K — » K + K Q \. 
Below we proceed by expanding f^, using the results 
in Eqs. (|3"5)) and (|3"9")) . in a power series in q for differ- 
ent limits in order to obtain a c \ and K e \ (note, however, 
in the expansion for the dielectric limit, for qL 3> 1, we 
also find terms odd in q). Note that since the tension 
and bending rigidities are identified through terms in the 
restoring force expansion which are proportional to h, 
second and higher order terms in h (see discussion at the 
end of the previous section) do not contribute to <j e \ and 
K cl . 

1. In the "small" screening length limit , n\L, k 3 L 3> 
1, a straightforward but lengthy expansion of 
in a power series in q assuming that the wave- 
length of the perturbation (= 2-K/q) is larger than 
the membrane thickness and the Debye screening 
lengths, qd,q/Ki,q/n 3 <C 1, gives 



= - [vciq 2 + K cl q 4 + 0(q 6 )] h. 



(41) 



B. Contribution to the free energy of the 
membrane 

Let us now investigate the electrostatic contribution to 
the membrane free energy. We first note that if the flu- 
ids surrounding the membrane are incompressible, then 
the pressures p° [occurring in the zero order terms in 



Eqs. (I38p and (|3"9")l ] adjust such that there is no net 
force (and hence no net movement of the membrane) in 
the z-direction; we will here consider such incompress- 
ible fluids and also assume the membrane to be incom- 
pressible. Nevertheless the investigation of the net force 
/(°) = /1I2 + /2-3 °n a membrane provides some insights 
into the electrostatic problem under consideration here, 
and is given in appendix 151 

Let us now proceed by considering the z-component of 
first order correction to the forces, = f^2 2 + /2-3> 
from p 1 ' one can obtain the work on the membrane un- 
der an undulation deformation of the shape, and thereby 



The explicit expression for the electrostatic contri- 
bution to the tension is 



£pAjfr| h+h + Ad/e 2 
2 (h +k + 2d/e 2 f 



(42) 



where, as before, the "rescaled" Debye screening 
lengths are l\ = (ei^i) -1 and l 3 = (£3/^3) -1 . We 
also introduced 



A<j) h+k + 2d/e 2 
2 Ji + l 3 + d/e 2 



(43) 



being the potential difference between the main 
parts of the two bulk fluids. Notice that a c \ gives 
a negative contribution to the tension. The fact 
that cr e i < originates from the fact that that 
the applied potential creates a net charge density 
on either side of the membrane surfaces, see Fig. 
[5] since ions of equal charge repel each other the 
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system would, for a compressible membrane, be 
able to decrease the free energy by separating the 
charges through a stretching of the membrane (i.e., 
by increasing the membrane area). For an incom- 
pressible membrane (as assumed here) the mem- 
brane is likely to respond to the electrostatically 
induced negative tension by an opposite increase in 
the membrane elastic contribution to the tension. If 
the magnitude of <7 e i exceeds the membrane elastic 
strength (tensile strength) a stretching instability 
occur. For the symmetric case (k = Ki = k 3 and 
e = Ei = £3) Eq. (|4*2")l becomes: 



(T e i 



o l + g m (Kd)-V2 

^(i+M/^)- 1 ) 2 

£ O £ 2 (A0 m ) 2 

2d 



(44) 



where we introduced the ratio e m = £2/2 between 
the membrane and surrounding medium dielectric 
constant (typically e m k 1/40, see [H, 0])- We 
notice that when s — £ m (Kd) _1 <C 1, i.e. for 
an effectively small screening length compared to 
the membrane width, the tension approaches er^; 
the dimensionless parameter s is commonly appear- 
ing in membrane electromechanical problems, see 
Refs. 0, 0j- From Eq. |g2J) we notice that for 
the asymmetric case we similarly have a c \ « crj?j for 
s-y = (£2/e 7 )(K(i) _1 <c 1, where 7 = 1,3. Since 
CTgj only depend on the membrane dielectric con- 
stant, membrane width, and the applied potential, 
<7 e ] is for large salt concentration independent on 
the properties of the surrounding medium (i.e. in- 
dependent on ki, K3, £1 and £3). The origin of this 
result is discussed below. 

The electrostatic contribution to the bending rigid- 
ity is 



K el = £ O A0; 



with coefficients 



b Q + bid + b 2 d 2 + 6 3 d 3 



M 4 



(h + h + 2d/£ 2 ) 3 



(45) 



bo = 

h = 

b 2 = 

h = 

b A = 



\(-(h^ 2 + hn~ 2 ) 

4 \£2 
6 £2 

4 1 
3iT 



We point out that K c \ > 0, i.e., the applied po- 
tential tends to make the membrane more rigid to- 
wards bending. During a bending deformation the 
induced charge density on one side of the mem- 
brane gets compressed, whereas the charge density 
on the opposite side gets expanded. The free energy 
changes of compression and expansion has different 
signs, but are in general of different magnitude. It 
is only for the case that all Debye screening charges 
are collapsed onto the surfaces (n — > 00) and zero 
membrane thickness d — > that the expansion and 
compression free energies are identical and K c \ = 
[see Eq. (|451) ]. Thus, loosely speaking, the smaller 
the "effective" membrane thickness (the membrane 
thickness including the Debye screening layer thick- 
nesses) the smaller is the bending rigidity. For the 
symmetric case (ft = Ki = K3 and e = £1 = £3) we 
write Eq. (|4"5)) according to 

K cl = ^(l + e^/sd)- 1 )- 3 
x (l + Aem(Kd)- 1 + 3£ m (l + £ m )(KcfT 2 

+e m (9/8 + 3£ m )M~ 3 + (9/8)£ 2 „ M" 



Kl = -£ o £ 2 d(A0 n 





(47) 



(46) 



We note that in the limit of small relative mem- 
brane dielectric constant e m — » as well as for 
small screening length compared to the membrane 
thickness, (red) -1 — > 0, we have K \ — > K®,. No- 
tice that, in practice, e m is always small w 1/40, 
and therefore we have K c \ — » also for "not 
too small" values for (kg?)" 1 . Similarly, we note 
from Eq. ([43]) that for the asymmetric case we 
have K c \ — > in the £2/67 — > and in the 
(Kjd)^ 1 — > (7 = 1,3) limits. For the above con- 
sidered limits the major part of the potential drop 
is across the membrane (since £2/£ 7 or (kg?) -1 is 
small), and hence the electric field is essentially 
zero everywhere except for the membrane region. 
Therefore, the membrane parameters play the dom- 
inant role in the expression for K c \ and <r c \ (notice 
that and er°j depend only on £o£ 2 , d and A<p m ), 
provided that the interactions of the induced sur- 
face charges do not occur through the surround- 
ing medium (again, certified if e 2 /£-y or (kg?) -1 is 
small) . Given that K e \ and <r c \ can only depend on 
£o£2 , d and A</> TO in the limits considered above one 
may use dimensional arguments to argue that the 
limiting results, and er^, for the bending rigid- 
ity and tension must (up to a constant prefactor) 
take the forms given in Eqs. (|33| and (|4"7|) . 

Fig. [3] illustrates the electrostatic contribution to 
the tension and membrane bending rigidity and its 
dependence on salt concentration (Debye screen- 
ing) for the symmetric case for simplicity. We see 
that the absolute value of electrostatic contribution 
to the membrane tension increases with increasing 
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salt concentration, i.e. for increasing k. We at- 
tribute this to an increase of screening charges (in 
a layer of decreased thickness) next to the mem- 
brane; the increased amount of charges will result 
in larger electrostatic repulsion between ions in the 
screening clouds [see discussion following Eq. (|43|) ]. 
In contrast to the effect on tension the electrostatic 
contribution to the bending rigidity decreases with 
increasing salt concentration (with K e \ approach- 
ing as k — ► oo). The reason behind this is that 
for bending properties the thickness of the Debye 
screening layers plays a role - a larger "effective" 
membrane thickness gives a higher bending rigidity 
[see discussion following Eq. (f45|) . and the sponta- 
neous curvature calculation in Appendix [C] . If we 
choose the potential difference between the mem- 
brane sides A<fi m = 100 mV, £2 = 2 and d — 2.5 nm 
we find that = 0.018fcsT (for room tempera- 
ture, k B T = 4 • 10~ 21 ). From Fig. [3] we see that 
K e i/K^ can become quite large for small k and, 
therefore, for sufficiently small salt concentration 
the bending rigidity K c \ can exceed the thermal 
energy fc^T; we therefore expect that the increase 
of the bending rigidity in the presence of an applied 
electrostatic potential predicted in this study can 
indeed be experimentally observed for sufficiently 
large A</> m and small salt concentrations (note how- 
ever, the below restriction on salt concentration 
due to assumptions in the Debye-Hiickel approxi- 
mation). 

Let us compare the results above for the symmet- 
ric case to the results obtained in Q . We find that 
the result for <7 e i given in Eq. (|4"4"|) agrees with the 
finite bilayer thickness, non-conductive membrane 
tension (there denoted by £j n + £ ut) obtained in 
Q (note that in Q the membrane thickness is de- 
noted by d, whereas we denote by d the size of a 
lipid monolayer, so that in our case the membrane 
thickness is 2d. Also note that, due to different 
boundary conditions at the electrodes, the V in Q 
should be equated with our A<p m ). Concerning the 
bending rigidity result, we note that no explicit ex- 
pression for K c \ was given in [?|, only a numerical 
value for a specific set of parameter values. We 
choose the same parameter values (A(/> m = 50 mV, 
d = 2.5 nm, e rn = 1/40, and 2dn = 7.4), but note 
that in order to get the expression for K e \ one must 
also choose the actual value of £2 (not just the ra- 
tio Eije), which was not specified in 0]. We choose 
the standard value £2 = 2 [3, 0| and then find that 
K ei = 0.00467 k B T, which is a bit less than half 
the value found in Ref. [7j. Since no explicit ex- 
pression for K c \ was given in Q it is difficult to 
comment on the nature of this discrepancy. Finally, 
using the approximate expression K c \ « for the 
parameters above we find that this approximation 
underestimates K c \ by merely 1 %. 

Here, a few words on the validity of the Debye- 



Hiickel approximation, used throughout this study, 
are in place. This approximation should work when 
the quantity / = j3q l (^ 7 (x) — 7 ), where 7 = 1, 3, is 
very small, i.e. / < 1 (see Sec. HI]) , but in practice 
the Debye-Hiickel approximation works well when- 
ever / < 1, see Ref jl9| . The maximum of / occurs 
at the membrane surfaces, see Eqs. (|15[) . (fl~6| and 
Fig. [21 and we find that for the k 7 L ^> 1 limit con- 
sidered here we have / = (3q l A(j}l 1 /[2(li+l 3 +d/e2)}- 
Using d » 2.5 nm, £2 ~ 2 and £ 7 ~ 80, the denomi- 
nator in the expression for I above is dominated by 
the djs2 term whenever k" 1 < 50 nm (in this limit 
also A<fi m Acj) m ); for such ion concentrations we 
have the following criterion for the validity of the 
Debye-Hiickel approximation 

/ = i^A^ < 1. (48) 

again involving the parameter s 7 = 
(£ 2 /£ 7 )(K 7 d) _1 . For a large potential differ- 
ence A<fi m — 100 mV we find that / < 1 when 
n^ 1 < 50 nm (using q 1 = 1.6 x 10~ 19 C and 
assuming room temperature). This means that 
the Dcbyc-Huckel approximation, which was made 
m Sec. HJ in the main text, works surprisingly 
well in general; the reason for this is the small 
value of £2/ £ 7 (« 1/40) guaranteeing that the 
major part of the potential drop occurs across 
the membrane and that, therefore, the potential 
drop across the electrolytes, for which we applied 
the Debye-Hiickel approximation, is modest, see 
Fig. [2 We point out that the dielectric limit, 
c buik,7 — * 0, does not rely on a Debye-Hiickel 
approximation and results to be given below 
are therefore valid for any value of the applied 
potential. We have shown above that our results 
for the small screening length limit (with the above 
explicit restriction) and in the dielectric limit are 
usually valid. We note, however, that there is 
in general an intermediate salt regime where the 
Debye-Hiickel approximation breaks down (for 
sufficiently large applied potentials). We leave the 
investigation of this intermediate regime for future 
studies. 

We finally point out that there are alternative ways 
of computing the membrane mechanical parame- 
ters. For instance one can calculate the tension by 
finding the integral of the deviation of the pressure 
profile from the value of the pressure far from the 
membrane. This approach is demonstrated in Ap- 
pendix giving the same result as Eq. (|4"!?]) for 
the tension. In that appendix the same type of ap- 
proach is also used in order to obtain the electro- 
static contribution to the membrane spontaneous 
curvature, see Eq. (|C7|) . 

2. in the dielectric limit, ki,k 3 — > 0, one can again 
perform a power series expansion in q using Eqs. 
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FIG. 3: The electrostatic contributions o c \ and K c \ to the 
tension and the bending rigidity and as a function of Debye 
screening in the small screening length limit (symmetric case), 
using Eq. (144[) and (|47[) . We express K c \ in terms of the (room 
temperature) thermal energy ksT and a e \ in terms of its infi- 
nite k value <t°j. The parameters used are listed in the figure. 
Notice that increasing salt concentration, i.e. increased k, 
leads to a decrease in the electrostatic contribution to the 
bending rigidity. In contrast, for increasing k the magnitude 
of the electrostatic contribution to the tension increases. 



(|38|) and (|39|) and assuming gd < 1. In addi- 
tion, there are two limits of interest depending on 
whether the wavelength perturbation is smaller or 
larger than L. 

For (i) 5L > 1 we find that 



/ 



(i) 



-[aq + <T cl q 2 + bq 3 + K cl q A + 0{q 5 )}h, (49) 



where 



-£o 



^r 1 



- 1 ) 3 



(50) 



and 



(Tel = -8(fe £2 



Acj) 



(51) 



and the higher order terms in q are more compli- 
cated functions of the dielectric constants and d. 
The term linear in q is negative and is expected 
to cause an instability for long wavelengths [20l |. 
We notice that the prefactor a [see Eq. (|50|) ] is 
proportional to the membrane polarization charge 
density (for a flat membrane) squared, a oc (p s ) 2 , 
see Eq. (|B4|) . The linear, non-analytic, q term 
may be interpreted as follows: for the asymmetric 
dielectric case the external potential induces an ef- 
fective net membrane polarization charge density 



p s . The membrane charges interact via the un- 
screened (there are no ions in the present limit) 
Coulomb interaction, giving rise to the (p s ) 2 pro- 
portionality for a; the non-analyticity of the free 
energy arises due to the long-range character of the 
Coulomb interaction (which decays as 1/r, where r 
is the distance between charges) . We point out that 
the linear q term does not depend on the membrane 
parameters, d and £2, which means that this insta- 
bility should exist for any interface between coex- 
isting fluids (fluid interface instabilities of a rather 
similar nature has been investigated previously, see 
for instance [(|). For the symmetric case £1 = £3 
we see that a — 0, but the third order term above 
is still present, b ^ 0, in general (see discussion 
below). 

For the case (ii) qL <C 1, i.e. the wavelength 
perturbation is longer than the distance between 
the electrodes, the first order force takes the form 
= -[c + (T c iq 2 + K cl q 4 + 0{q 6 )]h. In this 
limit the interactions between the membrane po- 
larization charges become short-range due to effec- 
tively small (compared to the wavelength of the 
perturbation) distance of the membrane to the elec- 
trodes. In this limit we thus do not have any odd 
q terms the effect of the applied potential is simply 
to give an electrostatic contribution to the tension 
and bending rigidity. The explicit expressions be- 
come somewhat complicated, but can be produced 
straightforwardly by a small q expansion using a 
symbolic mathematical software like Mathematica 
or Maple together the expressions for /i_2 and /2-3 
given in Eqs. ([55)) and 



For the symmetric case (e = e± = £3, K = K\ = K3) 
the lowest order term in the q expansion is the ten- 
sion term (q 2 term), both in the small screening 
length limit and for the dielectric limit. For both 
these limits we have that the electrostatic contribu- 
tion to the tension a c \ is negative. Thus when \<J e i\ 
becomes sufficiently large, i.e. large applied poten- 
tial, a membrane stretching instability can occur. 
We also find that in the symmetric dielectric limit 
there is a q 3 term present (for qL ;§> 1), where ex- 
plicitly we find [see Eq. (|4"9"1) ] 



&=2£ ££ 2 (^) 2 (£- 1 



(52) 



We notice that b oc P 2 , where P z is the membrane 
polarization per unit area, see Eq. (|B8[) . For the 
symmetric dielectric case the effect of the applied 
potential is to polarize the membrane and the g 3 - 
term is expected to originate from induced, un- 
screened, dipole-dipole interactions (which decays 
as 1/r 3 , where r is the distance between the dipoles 
[l5l |) within the membrane. 
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VI. SUMMARY AND DISCUSSION 

We have in this paper derived expressions for the elec- 
trostatic contributions to biomembrane mechanical pa- 
rameters (such as tension and bending rigidity) in the 
presence of an static applied potential across a mem- 
brane. The membrane was assumed non-compressible, 
non-conductive (membrane region described by Laplace 
equation) and surrounded by electrolyte solutions (de- 
scribed by the Debye-Hiickel equation). By solving the 
equations for the electrostatic potential and using the 
stress-tensor the forces acting on the membrane were ob- 
tained, which in turn were used to obtain the free en- 
ergy and the electrostatic contribution to the membrane 
mechanical parameters as a function of the applied po- 
tential, the salt concentrations (entering through the De- 
bye screening lengths) and the dielectric constants of the 
membrane and the solvents. Results of particular inter- 
est, that are found in this study, include: for (1.) the 
small screening length limit, where the Debye screening 
length is smaller than the distance between the elec- 
trodes, the screening certifies that all electrostatic in- 
teractions are short-range, leading to a free energy ex- 
pansion of the form ~ cr c \q 2 + K e \q 4 + 0(q 6 ) (where q is 
the wavevectors) , the main effects of the applied poten- 
tial are to decrease the membrane tension and increase 
the bending rigidity; explicit expression are given in Eqs. 
(|42)> and (|45|) . Our expression for the tension for the 
symmetric case reproduces the result in Q. In Q it was 
also found that an applied electric field gives a negative 
contribution to the tension. However in that study the 
medium surrounding the membrane was characterized by 
conductivities rather than Debye screening lengths, and 
it is therefore difficult to directly compare our results to 
theirs. For sufficiently large applied potentials the mag- 
nitude of the electrostatic contribution to the tension will 
exceed the maximum tension the membrane can sustain, 
leading to a membrane stretching instability. Possibly, 
this instability can result in the formation of pores and 
flow of ions through the membrane (in fact, the mem- 
brane tension is one of the key parameters in the model- 
ing of membrane electroporation dynamics [2l|). For (2.) 
the dielectric limit, i.e. no salt (for small wavevectors q 
compared to the distance between the electrodes), when 
the dielectric constants on the two sides are different, the 
applied potential induces an effective (unscreened) mem- 
brane charge density, whose long-range interaction causes 
a membrane undulation instability if the dielectric con- 
stants of the two bulk fluids are different; this effect is 
characterized by a negative term linear in q in the free en- 
ergy expansion, see Eqs (|49f and (|50|) . i.e. this term is of 
lower order in q than the tension term. Previous similar 
results include: In [|| the interface between two immisci- 
ble fluids of different dielectric constants was found to be 
unstable in the presence of a perpendicular electric field. 
The case of stiff (charged) DNA with bound, but mobile, 
counter ions was investigated in 0] and a shape instabil- 
ity found (supposedly leading to a DNA condensation). 



If the two dielectrics on each side of the membrane are 
identical, we found that then the applied electric field 
will give a negative contribution to the tension. Hence, 
if the applied potential is sufficiently large a membrane 
instability occurs also for the (dielectric) symmetric case. 

We quantified the validity of the Debye-Hiickel approx- 
imation (used throughout this study) and showed that 
our results are in general valid in the small screening 
length limit as well as in the dielectric limit. However 
for "small" , but non-zero, salt concentration and large 
applied potential the Debye-Hiickel approximation is no 
longer valid and one needs to consider the full Poisson- 
Boltzmann equation. It remains a future challenge to 
solve the full Poisson-Boltzmann problem in order to find 
expressions for the free energy for arbitrary salt concen- 
tration, and, in particular, to investigate in more detail 
the nature of the onset of the predicted membrane insta- 
bility (via the negative term linear in q in the free energy) 
as salt concentration is lowered. 

Possible applications of the result for the small screen- 
ing length limit above would be to lipid membranes where 
a potential difference across the membrane is enforced 
by ion pumps incorporated in the membrane. Changes 
in membrane rigidity might then be observed in mi- 
cropipette or video microscopy experiments, if the screen- 
ing length and the membrane potential are large, see Fig. 
131 Another possible experiment would be to observe the 
structural change of domains in a multi-component mem- 
brane using fluorescence correlation spectroscopy (FCS) 
as the membrane potential in a patch clamp experiment 
is altered; one might be able to observe a change from a 
phase with caps to one with stripes or buds as the effec- 
tive bending rigidity is changed by the applied potential 

M 

For the cases discussed above where a membrane in- 
stability occurs, the system is expected to be driven from 
the quasi-flat shape into a new equilibrium configuration. 
Our perturbation analysis cannot in general say anything 
about this new configuration. However, for the small 
screening length limit we above speculated that the neg- 
ative electrostatic contribution to the tension (for large 
applied potentials) could lead to electroporation and a 
corresponding flux of ions through the membrane. It 
may also be speculated (similar to the studies in Refs. 
[2| and [9() that the electrostatically induced new equi- 
librium configuration under certain conditions could be 
a spherical membrane (a vesicle); it is known that vesi- 
cles can be created in laboratories in a process known 
as electroformation [f| under the application of electric 
fields. We hope that the present theory will stimulate 
further work directed towards the controlled estimation 
of vesicle sizes as a function of electrostatic parameters, 
e.g. potential differences and electrolyte concentration. 
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APPENDIX A: DERIVATION OF BULK 
EQUATIONS 

In this appendix we show how the electrostatic equa- 
tions and the stress-tensor used in the main text can be 
derived as Euler-Lagrange equations from a free energy. 
Suppressing the region index 7 for convenience we can 
write the appropriate free energy as 



1 



G = J d 3 x 

+ ^c l (fc B T(ln 



^e e\V<Z>(x)\ 2 +/?$ -p 



CTotal 



(Al) 



Here p = p{x) is the local pressure, c % — c l (x) are the 
local concentrations of the ions, CTotal = CTotai(a^) is the 
total concentration of molecules of all kinds (including 
water) and // are (constant) chemical potentials for the 
ions. The logarithmic term corresponds to the entropy 
of mixing. Some of the Euler-Lagrange equations of this 
free energy can be found by demanding stationarity when 
varying the ion concentrations 



Sc l 



9q l + k B T\n- 



CTotal 



+ M 1 



Solving for the ion concentrations we find 



CTotal exp 



q l $ + ^} 



(A2) 



(A3) 



Inspection of the expression for the charge density p — 
5^ <?V with c 1 given by Eq. (|A3j) reveals that there is a 
unique value of $ for which p vanishes. We define <fp to 
be this "equilibrium value" , i.e. 

pU =0 o = 0. (A4) 

The ion concentrations at "equilibrium" is labeled by 

c bulk = C 1$ = 0O ■ (A5) 

The Euler-Lagrange equation for <E> is 







5G 



e eV 2 $ + p, 



(A6) 



which is simply the Poisson equation. Insertion of Eq. 
13[) and Eq. (|A5[) gives the Poisson-Boltzmann equation 



<f ($ - 4>°) 



= 0, (A7) 



which when linearized gives the Debye-Huckel equation, 
Eq. © in the main text. 

The force-balance in the bulk regions can be found by 
demanding that G should be stationary with respect to 
moving all fluid elements, particles and fields by an in- 
finitesimal position dependent distance 5x — 8x(x). De- 
noting fields after the move by a prime one has the new 
fields 



) = (A8) 
p'(x') = p(x), (A9) 
c l \x') = (1- V-5x)c l {x), (A10) 

+ Sx. The free energy after the move is 



where x! = 
G' = 



d 3 x' 



- eo e|V'$'(x')| 2 



k B T In ■ 



/£>'$' - p' 



-Total 



and using 



d 6 x' = {l + V-8x)d c 



x. 



d 



E 



dSxj 
dxi 



d 

dxj 



one finds that the change in free energy is 

dSxj 



5G = G 



f -G=[d 3 x y^r,. 



dxi 



(All) 

(A12) 
(A13) 

(A14) 



where is the stress tensor for the system 



(89 89 1 ^ 

Tij - £0£ [dx~dx-~2 Sl ^ 



dx k dxk 



pSij. (A15) 



Note that the first two terms on the right hand side of 
Eq. (|A15[) are simply the classic Maxwell stress tensor. 
Performing a partial integration in Eq. (|A14[) one can 
see that the condition of stationarity of G implies the 
conservation of stress 



dxi 



0. 



(A16) 



Combing this with Eq. 
physically revealing form 



61) one arrives at the more 



(A17) 



pV$ - v> = 0, 



i.e. electric forces should be balanced by changes in hy- 
drostatic pressure. Finally, since we know the behavior of 
the charge density p as a function of <& from Eq. (|A3|1 we 
can integrate Eq. (|A17|) to find the pressure as a function 
of $ 



exp 



1 



(A18) 
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The quadratic version of Eq. (|A15|) with Eq. (|A18|) in- 
serted is Eq. (f3T>|) of the main text. Note that the pres- 
sure can also be written p = p° + ^T^fc 1 — c l hulk ), i.e. 
there is an osmotic contribution to the pressure when ions 
are "confined" by the electric potential (see for instance 

Hi). 



APPENDIX B: 



TOTAL NET FORCE ON A FLAT 
MEMBRANE 



In this appendix we investigate the total force acting 
on a flat membrane in some detail. 

From Eqs. (03), (03) and Eq. (TTTJ) we obtain the 
following explicit expression for the total force (per unit 
area) acting on a flat membrane: 



(0) 



f 



£o 
' 2 



AO) 

Jl-2 
Pi 



AO) 

J2-3 



-l e - Kl L 



-k 3 L 



( 1 + e - Kl L)2 + , 



-n 3 L\2 



r 2 (bi) 



and T = l/\g(K\)h + g(K 3 ) h +d/e 2 ], and h = (eiKj.) 1 
and I3 = (£3K3) _1 as before [the function g(q) is defined 
in Eq. (|19|) ]. If the fluids surrounding the membrane are 
incompressible, then the pressures p® adjust such that 
the above force vanishes. In fact, in this study we assume 
that both the membrane and the surrounding fluids are 
incompressible. Let us, however, in order to gain some 
physical insights, investigate different limits of the non- 
pressure part of the above force (i.e. we take p\ = p®). 

1. In the "small" screening length limit, K\L, k 3 L 3> 
1, the non-pressure part of Eq. (|Blj) becomes: 



/<°> « 0. 



(B2) 



The force is zero whenever the screening length is 
small compared to the distance between the elec- 
trodes. Thus it is not possible to get a net force 
on the membrane solely by having different con- 
centration of ions (free charges) on the two sides. 
That there is no net force on the membrane can be 
related to the fact that their is no net free charge 
around the membrane, which can be understood 
from Gauss' law, J D ■ ndS — Q[ rcc , where D is 
the displacement field and Qf r0 c is the enclosed 
free charge: since far from the membrane on both 
sides the electric field is zero (in the here consid- 
ered limit, see Eqs. (fT5j). (flUI) and Fig. (3) so is 
the displacement field, and applying Gauss law (for 
a large "pillbox" enclosing the membrane and the 
Debye screening layers) one finds that <5f re c = 0, 
as well as no net charge due to changes in polar- 
ization of the dielectric media. The effective net 
charge for the membrane and the Debye screening 
layers is hence zero and there will be no net force 
due to these charges. 



2. In the dielectric limit, K\, K3 — > 0, we find 



^ = -2e (#f 



2L 



2(d/L)e 



(B3) 



for the non-pressure part of the force. When there 
are no free charges, as considered here, the Gauss 
law argument above does not apply. One then 
needs to also take into account the bound charges, 
which create a polarization surface charge den- 
sity pf_ 2 = — (P 2 - Pi) • z (see [H[, chap. 4) 
on the region 1-2 interface. Similarly, for the re- 
gion 2-3 interface we have a polarization surface 
charge density p|_ 3 = (P 2 — P3) ■ z. From the 
fact that P = D — EqE and that the normal com- 
ponent of the -D-field is continuous we find that 
P1-2 =sq(E 2 -Ex)-z and p^_ 3 = -e (E 2 -E 3 )-z, 
i.e. the polarization surface charge density is de- 
termined by the jump in the electric field at the 
interface. Using the explicit expressions for the po- 
tential $(2) given in Sec. IHII and that z-component 
of the electric field is E z = —d<fr/dz we find 



P1-2 
P2-3 



~£0-Ez,appl( e l 1 
£0-Kz,appl(£ 



^/(el 1 + e^ 1 ) and 



-3 



£ 2 1 )/( £ 1 1 + e 3 *)' with tne 

applied field being E z>app \ — Acj>/ (2L) where we as- 
sumed a small membrane thickness d/L <C 1. The 
effective membrane charge density p s = /Oi_2+/°2-3 
then becomes proportional to the dielectric asym- 
metry between region 1 and region 3, explicitly 



(B4) 



£0-^js,appl _i 



where the proportionality p s cx £3 — follows 
directly from the fact that the jump in the electric 
field is proportional to ea/ei (£2/ £3) for the region 
1-2 (region 2-3) interface (see discussion at the end 
of Sec. IIIII and Fig. (3)- Thus, whenever e\ ^ £3 we 
have that the membrane has an effective net charge 
of induced bound charges, which when put into the 
applied electric field E z ^ vv \ gives rise to the force 
above; note that the force is proportional to p s as it 
should (for d/L < 1). We also note that Eq. (fB3|) 
for d — > can be written 

/W = ^|-lWi(*-£i) 2 , (B5) 

which can be compared with for instance the for- 
mula in [l8| for the pressure decrease in a fluid at 
the fluid-air interface due to a normal electric field 



z-Ei = 
on the fluid side. 



Ac/) ^i 1 



(B6) 
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3. For the symmetric case (ei = £3 = e, n\ = K3 and 
Pi = P3) the non-pressure part of Eq. (|B 1 1) is zero, 



/ (0) - 0, 



(B7) 



f(0) 



on the two region boundaries should vanish, /j_ 2 



and / 2 (0) 3 

From either Eq. (J38J) or Eq. j39|) one finds 



(or equivalently: Tz° z \z) — — p° for all z) 



since for the symmetric case there is no net surface 
charge, and hence no net force on the flat mem- 
brane. We point out, however, that the membrane 
gets polarized. In particular, the polarization per 
unit area in the z-direction (using the results above) 
is 



— eosE z . app i(e 1 — e 2 1 )d 
in the dielectric limit. 



(B8) 



APPENDIX C: ALTERNATIVE APPROACH TO 
TENSION AND SPONTANEOUS CURVATURE 

In this appendix we will demonstrate an alternative ap- 
proach to deriving the tension in the membrane, which 
can also be used to obtain the spontaneous curvature 
induced by the electric field. This approach consists of 
calculating moments of the deviation of the pressure pro- 
file from the value of the pressure far from the membrane. 
We will therefore in this appendix only study the small 
screening length limit where L — ► 00, such that the pres- 
sure approaches a well defined pressure away from the 
membrane. 

The tension is obtained as the integral of the lateral 
pressure profile deviation, or equivalently the excess lat- 
eral stress, of the planar membrane 241] . Choosing the 



diagonal x-component of the stress tensor to represent 
the lateral stress one has the precise formula [27 



dz 



T (o) 



(CI) 



where we have used that the pressure on the two sides 
of the membrane should be identical p° — p\ — p®, since 
the zeroth order force in the small screening length limit, 
Eq (|B2|) . vanishes. In terms of the zeroth order solution 
given in Eqs. H3), (fT5")l and (fit?]) we obtain 



so 



e 1 K 1 A 1 - — e 3 K3^ 3 ~ e Q e 2 dA 2 



-2d(p° 2 -p°) 



(C2) 



To obtain the final result we need to find p 2 . To do this 
we note that in a flat equilibrium configuration the force 



This gives 



P°2=P°+ 2 £ o £ 2 A l 



60 - - <» £o -e 3 n 3 Al 



2e e 2 dA? 2 



(C3) 



(C4) 



This expression is identical to the tension given by Eq. 
(|4"2")l in the main text. 

An advantage of the approach of integrating the stress 
profile is that we can obtain the change in spontaneous 
curvature Co induced by the electric potential. If we in- 
clude a spontaneous curvature in the Helfrich free energy 
G given at the beginning of Sec. IV Bi writing it as 



G= dA 



-K{2H) 2 - KC Q 2H 



(C5) 



and, as before, calculate the force / rs = —SG/5h, then 
the spontaneous curvature will drop out at linear order in 
h and we would end up with our previous expression for 
the force where Co is not present; thus it is not possible 
to obtain the spontaneous curvature using the approach 
of the main text. However, just like for the tension, the 
spontaneous curvature can be obtained from the lateral 
stress profile, namely as the negative first moment of the 
lateral stress profile, sometimes also called the bending 
moment. The formula is ]2M 



KC 



dz z 



Ti°J - (-/) 



(C6) 



and insertion of the zeroth order solution gives 



KCo = -^A\{1 - 2 Kl d) - ^f-A\{l - 2K 3 d) 



or explicitly 



KCo = j (4^) (Wr(l ~ 2« x d) - £3/1(1 ^ 2 K3 d))r 2 , 

(C7) 

where T = l/[h + I3 + d/e 2 ], h — (eiKi) 1 and ^3 = 
(e3«3) _1 - Note that Co vanishes in the symmetric case 
where n\ — K3 and e\ = e 3 as it should. 
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